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We study superfluidity of supersolid phases of dipolar Bose gases in two-dimensional optical 
lattices. We perform linear stability analyses for the corresponding dipolar Bose-Hubbard model in 
the hardcore boson limit to show that a supersolid can have stable superflow until the flow velocity 
reaches a certain critical value. The critical velocity for the supersolid is found to be signiflcantly 
smaller than that for a conventional superfluid phase. We propose that the critical velocity can be 
used as a signature to identify the superfluidity of the supersolid phase in experiment. 
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I. INTRODUCTION 

New possibilities to explore exotic quantum phases 
have been pioneered by recent experimental advances 
in creating dipolar ultracold gases [ll, [H, such as the 
reahzation of a condensate of ^^Cr atoms with strong 
magnetic dipole moments Q and heteronuclear po- 
lar molecules [H, Thanks to the long-range nature 
and anisotropy of the dipole-dipole interactions, various 
quantum phases have been predicted to emerge, includ- 
ing fermionic superfluids (SF) with p-wave pairing [2 , [71 , 
Haldane-Bose insulators Q, and supersolids (SS) |9l-[l4j. 

Of particular interest are SS phases, in which both 
diagonal (crystalline) and off-diagonal (superfluid) long- 
range orders coexist [l^. Although non-classical ro- 
tational inertia, one of the signatures of su per fluidity, 
was experimentally observed in solid helium [l6[, it has 
been more reasonably interpreted by other mechanisms, 
such as superfluidity of grain boundaries and disloca- 
tions [TtI . fisj . On the other hand, quantum Monte 
Carlo (QMC) simulations have shown the presence of SS 
phases in Bose-Hubbard systems with long-range inter- 
actions [ll,[ll,[il[23l. Since QMC analyses of the Bose- 
Hubbard model quantitatively agree with experiments of 
ultracold gases in optical lattices, the SS phases are ex- 
pected to be found in the context of dipolar Bose gases 
loaded into optical lattices. 

In order to verify the existence of SS in experiments 
of ultracold gases, one has to clarify observables to iden- 
tify the superfluidity and the crystalline order of the SS 
phases. It is well-known that the crystalline order can be 
identified by the static structure factor, which has been 
observed in cold atom experiments using the Bragg scat- 
tering techniques [U . On the other hand, a sharp 
interference peak in the time-of-flight image following the 
expansion of a gas is often used as an indirect indication 
of the superfluidity of Bose gases in optical lattices [l^ . 
However, the sharp interference peak identifies the pres- 
ence of a Bose-Einstein condensate, but does not nec- 
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essarily mean the superfluidity. For instance, although 
a non-interacting Bose gas forms a condensate at suffi- 
ciently low temperature, it is not a SF in the sense that 
its critical velocity is zero [13]. Moreover, while most 
previous theoretical work calculated the SF fraction as 
a characteristic of the superfluidity [l^, [3 [H, US], so 
far no experiment has succeeded in measuring the SF 
fraction in cold atom systems [25j . in contrast to helium 
systems where the SF fraction, corresponding to the non- 
classical rotational inertia, can be easily measured with 
a torsional oscillator [l^, [l^ ■ Instead, the superflui dity 
of weakly- [g^l and strongly- interacting Bose gases [28j |. 
and fermionic SF across the BEC-BCS crossover [2^ has 
been demonstrated in a moving optical lattice by mea- 
suring the critical velocity above which superflow breaks 
down. 



In this paper, we propose that the superfluidity of 
SS phases can also be experimentally identifled using a 
moving optical lattice. Performing linear stability anal- 
yses for polarized dipolar hardcore bosons in a two- 
dimensional (2D) moving optical lattice, we show that 
superflow of SS states is stable until the flow momen- 
tum exceeds a certain finite value. It is found that the 
critical momenta for the dynamical instability in the SS 
phases are remarkably smaller than that for a standard 
SF with no density wave order. We argue that the critical 
momenta can be experimentally measured with currently 
available techniques. 



The remainder of the paper is organized as follows. In 
Sec. ini we introduce our model Hamiltonian describing 
hardcore bosons with dipole-dipole interaction in a 2D 
optical lattice. In Sec. IIIIl we explain our formulation of 
the problem based on a mean-field theory. In Sec. IIVI 
we calculate the ground-state phase diagram. In Sec. |Vl 
we perform linear stability analyses to obtain the critical 
velocities for Landau and dynamical instabilities in the 
SF and SS phases. In Sec. IVIl we summarize our results 
and briefly discuss two recent experiments exploring SS 
phases in cold atom systems [s^ [3l| . 
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II. MODEL 

We consider a system of N bosons interacting with on- 
site and dipole-dipole interactions in a 2D optical lattice. 
The dipoles are assumed to be polarized to the direction 
perpendicular to the lattice plane. This system can be 
well-described by the dipolar Bose- Hubbard model 



H = -jJ2ia]ai+h.c.) + -J2njinj -I) 
U-i) 

3<l 3 



(1) 



where is the boson creation operator at site j, hj = 

a^Oj, J is the hopping, and U is the onsite interaction. 
The chemical potential /i controls the filling factor n = 
N/M, where M is the total number of lattice sites, (j, /) 
represents nearest-neighbor pairs of lattice sites. The 
long-range part of the dipole-dipole interaction is well- 
approximated as Vji = Vd'^lvj — ri\~^, where jx and jy 
are integers and d is the lattice spacing. In experiments, 
the ratios J/V and U/V can be controlled by varying the 
lattice depth and using the Feshbach resonance [J| . 

In the hardcore boson limit (U — oo), Eq. ([T]) can be 
mapped onto the following spin-1/2 Hamiltonian, 

= -jJ2iS+S[- + h.c.) + J2 y^S-S^ - /i^ 5|,(2) 
0,0 i<i i 

where the spin operators are related to the operators of 
the hardcore boson as 5*" = aj and = -hj — 1/2. 
h = ji — 2C+V is the shifted chemical potential, where 



(±1)^ 



These coefficients include the long-range nature of the 
dipole-dipole interaction and their numerical values are 
C+ ~ 2.258 and C_ ~ 0.6615. The dynamics of the spin 
model Eq. ([2]) is described by the Heisenberg equation of 
motion for Sj' (we set h — 1), 



hS' 



(4) 



where {I) represents the nearest-neighboring sites to site 
j- 

III. MEAN-FIELD THEORY 

Our formulation of the problem is based on a mean- 
field theory, in which a many-body wave function is ap- 
proximated as a product of the local spin coherent states. 



1$ 



mf; 



n 



COS^|t)j 



e 2 sin^lDj 



(5) 



where 6j and tpj are the elevation and azimuthal angles 
of the spin direction at site j . Replacing 5'| and with 
(S";) = IcosOj and (5+) = \e'^^ smOj in Eqs. @ and 
we obtain the mean-field energy "Ho = (-^^s) given by 



cos 9i cos 9, 



jl cos Oj 



(6) 



and the classical equations of motion for 9j and ipj : 



^ = J ^ sin 6*; sin ipji , (7) 



(0 



d(pj J sin 9i cos 9j 
dt ^ sin 9 a 



cos Lfji 



l^Vji cos 9i 



h,{8) 



where ipji = tpj — ipi. Equation ([7]) is the continuity 
equation while Eq. ([5]) corresponds to the Josephson's 
acceleration equation. The hardcore boson density rij = 



the density 



|(aj)| , and the phase 



arg((aj)) of the condensate wave function, and the spin 
angles are related by rij — {cos9j + l)/2, n^™ = rij(l — 
rij), and = —^j- 

Considering small-amplitude oscillations around a 
steady solution, we write the solution of Eqs. ([7]) and 
dH) in the form 



9,{t) = 9, + 89,e-^-\ ^,{t) = + S^,e 



(9) 



where uj is the frequency of the normal mode. Substi- 
tuting Eq. (|9]) into Eqs. ([7]) and dH) and neglecting the 
terms higher than the first order with respect to S9j and 
6(pj, we obtain the equations for a steady state 



sin 9i sin (pji 

(0 



0, 



(10) 



■^-^ sin 9i cos 9j 



sm9i 



',1 COS0z,(ll) 



and the linearized equations for fluctuations 
— iuj59j = [^^ji sin 9i cos ipji+69i cos 9i sin (pji^ ,(12) 



(0 
(0 
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cos t/ , COS ( 
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„ sin 9i 
sm ( 



-^^,,(50, sin ^,.(13) 

1^3 



The excitation energy uj calculated from Eqs. p2|) and 
(|13p coincides with that calculated by linear spin-wave 
theory [s^] . Stability of a steady solution can be discrim- 
inated by ui. The appearance of excitations with a; < 
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FIG. 1: (color online) (a) Ground-state phase diagram of the 
hardcore Bose-Hubbard model with the dipole-dipole inter- 
action in the (1// J, /i/ J)-plane, where MI, SF, CSS, CS, and 
SS2 phases are present. The dashed red line represents the 
contour of n = 0.4. (b) The filling factor n as a function 
of h/ J . The dotted red lines locate the boundaries between 
the different phases. The dipolar interaction is fixed to be 
V = 3. 5 J, which is indicated by the dotted blue line in (a). 
(I) and (II) Schematic pictures of the CSS and SS2 phases. 



signals the Landau instability (LI) , while the appearance 
of excitations with Im[a;] ^ signals the dynamical in- 
stability (DI), which means exponential growth of the 
fluctuations in time. The linear stability analyses on the 
basis of Eqs. ([T^ and ([T51) allow us to calculate the crit- 
ical velocity of superflow. 

Previous theoretical analyses have shown that the 
mean-field theory fails to describe even qualitatively 
the ground-state phase diagram of the hardcore Bose- 
Hubbard model with nearest- and next- nearest- neighbor 
interactions due to strong quantum fluctuations [20l Is^ . 
More specifically, while the mean-field theory predicts 
the presence of a stable checkerboard supersolid (CSS) 
phase [131, this SS phase is unstable towards phase sepa- 
ration according to accurate QMC simulations. However, 
the mean-field theory is qualitatively valid for our dipo- 
lar hardcore Bose-Hubbard model of Eq. ([2]), because 
the long-range nature of the dipolar interaction signifi- 
cantly suppresses quantum fluctuations [s^ . We will in- 
deed show in the next section that the mean-field phase 
diagram for Eq. ([2]) qualitatively agrees with the recent 
QMC results of Ref. • We do not argue that the mean- 
field theory can provide quantitatively correct results, but 
that it is useful for gaining qualitative features and ana- 
lytical insights of the critical velocity. Notice that so far 
the QMC methods have not succeeded in calculating the 
critical velocity. 



IV. PHASE DIAGRAM FOR K = Q 

Within the mean-field theory, let us first calculate the 
ground-state phase diagram of Eq. ^ in the case that 
the system does not have a supercurrent, i.e. (pj is 
constant. Without loss of generality, we can consider 
the spins ordered in the XZ plane {(pj =0). To ob- 
tain analytical expressions of phase boundaries, we as- 



sume that the stationar 
jx + jy € even(odd) [3^ 



solution satisfies 0, 



for 



34| . Under this two-sublattice 
ansatz, we can describe checkerboard solid (CS), SF, 
CSS, and Mott insulator (MI) phases. The CS phase, 
termed as the Neel state in the language of the spin 
model, is an incompressible insulating phase at half fill- 
ing with a density wave order whose ordering vector is 
kjr = {TT/d,TT/d). This phase is favored in the region 
of J, \h\ ^ V, where the antiferromagnetic Ising term 
is dominant in Eq. ([2]). The SF phase is characterized 
by uniform density and finite condensate fraction, i.e. 
rij = n and n^°" 0. The latter condition reflects 
the existence of the off-diagonal long-range order. This 
phase corresponds to a canted ferromagnetic state in the 
spin system and is favored when J/V is large or 
is moderately large. The CSS phase possesses both the 
SF and checkerboard density- wave orders. Recent quan- 
tum Monte-Carlo simulations of Eq. ([2]) in Ref. [HI have 
shown that the CSS phase is indeed present in the inter- 
mediate region between the SF and CS phases. MI is an 
incompressible phase with n = or 1 , which corresponds 
in the spin language to a fully polarized magnetic phase 
in a strong magnetic field \h\. In terms of the spin an- 
gle, the conditions for the different phases to emerge is 
as follows: 

cos 6a ^ — cosOb = 1, CS, 

eA = OB and sin 9 A ^ 0, SF, 

Oa ^ Ob and smOA ^ and sinews =^ 0, CSS, 

cosOa = cos 6*5 = ±1, MI. (14) 

Minimizing the mean-field energy. 
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.2 , 



— JsmQA sin^^B + -r(C+ + C_) cos^a cos9b 



+ ^{C+~ C_)(cos2 6 A + cos^ Ob 
8 

-^{cosOa + cosOb), 



(15) 



with respect to 9a and 9b ^ we obtain the phase dia- 
gram in the (V/ J,h/ J)-p\a.nc as shown in Fig.[TJa). Re- 
flecting the particle-hole symmetry of the hardcore Bose- 
Hubbard model, the phase diagram is symmetric with 
respect to the line h = 0. 

In Fig. [Ija), it is seen that the CSS and CS phases 
are present when V/J> 2/C_. Increasing h/J from the 
n = MI region with a fixed value oi V/J > 2/C_, 
the system exhibits the continuous transitions to SF at 
h = /igp, CSS at h = /less, ^^'^ at h = in order. 
This behavior is clearly illustrated in Fig. [Ub), where 
the filling factor n is plotted as a function of h/J. The 
critical values of h for these transitions are given by 
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±iAJ + 2C+V), 
±2(C,V^ + 2J)^^_^^^^, 
±2V(C_y)2 -4J2. 



(16) 
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FIG. 2: (color online) Schematic pictures of the SF and CSS 
states with superflow in the spin representation. The arrows 
represent the local direction of the spins. The radius of the 
circles denotes the projection of the spins onto the xy-plane, 
which corresponds to the local condensate density n^"". The 
red dashed lines denote the projection onto the z-axis corre- 
sponding to the local density Uj. 
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FIG. 3: (color online) Stability phase diagrams of hardcore 
bosons flowing in a 2D lattice with quasi-momentum K = 
(iv, 0). In (a) the regions of stable SF, stable CSS, stable 
SS2, Landau instability (LI), and dynamical instability (DI) 
are located when V/J is varied for n = 0.4, while in (b) those 
are located when the filling factor n is varied for V — 3. 5 J. 
The dotted blue, solid red, and dashed green lines represent 
the critical quasi-momenta for LI, DI caused by phonons, and 
DI caused by roton-like excitations. 



We note that replacing C±V with 2{Vi ± V2), Eq. ^ 
coincides with the critical values of h obtained in pre- 
vious work for the hardcore Bose-Hubbard model with 
the nearest-neighbor interaction Vi and the next-nearest- 
neighbor interaction V2 [H, Ull . When V/ J is increased 
further, there emerge different solid and SS phases in ad- 
dition to the phases described above. For instance, allow- 
ing for the four-sublattice density modulation, we calcu- 



late the boundary to the SS2 phase as shown in Fig.[lja). 
The SS2 phase is sketched in Fig. [Tfll). We do not push 
our calculations into the region ofV > 5.2 J. where other 
solid and SS phases are present, because our purpose is 
to investigate superfluidity of the SF and CSS phases. 
Notice that there also exist numerous meta-stable states 
in the region of large V/J [IBl, which make experimental 
investigation of the ground-state phase diagram practi- 
cally very difficult. 

V. EXCITATION SPECTRA AND CRITICAL 
VELOCITY 

Having established the location of the SF and CSS 
phases in the phase diagram, we next study stability of 
superflow in these phases by means of a linear stability 
analysis. Let us consider that the optical lattice confining 
hardcore bosons is moving at a constant velocity v. In 
the coordinate system where the lattice is at rest, the SF 
component of the hardcore bosons is flowing with quasi- 
momentum K = — mv [13, [11], where m is the particle 
mass. 

In the SF phase, a current-carrying solution of 
Eqs. (fTO)l and (fTTj) is given by 6j = Oq and (pj = 
— K • Vj, where Oq is related to the filling factor as 
cos^o = 2n — 1. This state is sketched in Fig. [2][a). 
The current carried by this state is given by j = 2n(l — 
'^)JJ2m=x,y^m sm{Kmd), where {K^,Ky) = K, and 
and Gy represent the unit vectors in the x and y direc- 
tions. Inserting this solution into Eqs. ^ and (fTTj) . we 
obtain the energy per particle ck = (Ho + h^-{Sj))/N 
and the chemical potential /ik as functions of K: 



EK = -4(1 - n)J7K 



(1 - 2nf 
2n 



c+v, 



= 2(271 - 1)(2J7K + C+V). 



(17) 



(18) 



where 7k — J2m=x y^'^^i^md)/2. Notice that ek and 
/iK satisfy the thermodynamic relation, /ik — (^(^"^i^^ 



dN 



Solving Eqs. ([T^ and p^ . we obtain the excitation 
spectrum. 



WK(q) = (2-4n)J(7q_K-7q+K)+[4J(27K-7q+K-7q-K){2J7K-(l-2n)2j(7q+K+7q-K) + n(l-n)T^(q)}]^(19) 

I 



where q is the quasi-momentum of the excitation 
and V{q)/V = Y.'^=i {cos{aq^d) + cos{aqyd)} + 
cos{a^q,rd) cos{ayqyd). Stability of 

superflow in the SF state can be judged by Eq. 

We assume that the current is flowing in the x-direction, 

i.e. K = {K, 0), and depict the stability phase diagrams 



in the (n, Kd)- and {V/ J, Kd)-p\aiies in Fig. [31 Stabil- 
ity of the hardcore boson system deep in the SF region, 
V'<JorO<7i<C 0.5, is analogous to that of 
softcore boson systems described by the CP mean-field 
theory [s^, [IHl in the following way. With increasing K, 
excitations with CL'K(q) < appear at a certain value of 
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FIG. 4: (color online) Excitation spectra i^(q) in the CSS 
phase for V — 3.5J, n — 0.4, and different values of K in 
the a;-direction, where K — (a), K = -ft'css (b), and K — 
0.83/d > -K^css (c)- (d) and (e) are magnification of (b) and 
(c) focused on the region of \qx\ <^ 1/d and qy — 0, where 
excitations with negative and complex energies arise. In (f), 
the imaginary part of i^(q) corresponding to (e) is shown. 



the quasi- momentum, K = -R'gp, signaling LI. In Fig. |31 
Kgp is plotted by the dotted blue line separating the 
stable SF and LI regions. When K > K^p = n/{2d), 
long- wavelength phonons cause DI. This DI reflects the 
fact that the effective mass in the x-direction, defined by 
(to^)~^ = gj^^ , is negative, resulting in the imaginary 

sound speed in the a; direction c|-- = (kk^k)"^^^- Here, 
kk is the compressibility. 

In the SF region close to the boundary with the CSS 
phase, the excitation spectrum a;K(q) has a roton-like 
minimum at q = k,n- [32] ■ When K increases in this re- 
gion, the roton-like excitations cause DI, which signals 
the transition to the CSS phase [39l - l4^ . before LI oc- 
curs. The critical value of K for this DI is plotted by the 
dashed-dotted green line in Fig. [31 Previous theoretical 
work for the system of hardcore bosons with only nearest- 
neighbor interaction has predicted that superflow can be 
destabilized by the roton-like excitations, but that the 
resultin g C SS state with superflow is dynamically unsta- 
ble [iol. l4ll|. In contrast, we will show below that the 
flowing CSS state can be stable in our system of dipolar 
hardcore bosons. 

In the CSS phase, solving Eqs. ([T0|) and ([TT|) within 
the two-sublattice ansatz, we obtain a current-carrying 
solution, 

(pj ~ — K • Vj, cos 9 A = 4?! — 2 — cos 9b, 



cos6'b = 271 - 1 ± 2 Jn2 - n + i - ^- — (20) 
V 2 2pK 



where /3k = \/{C-Vy — (2J7k)^- A schematic picture 
of this state is depicted in Fig.[2jb). This state possesses 
the current j = 2|1 - 2n\J^jKPK J2m=x,y ^™ sm{Kmd). 
Substituting Eq. 1^ into Eqs. © and ([IJ), we obtain 



the energy per particle and the chemical potential. 



n 



K 



2n 



V, (21) 



/iK = sgn(2n - 1) X 2/Jk - 2(1 - 2n){C+ - C^)V. (22) 

In Fig. [H assuming K = {K, 0) again, we show the ex- 
citation spectra for V = 3.5J, n = 0.4, and the different 
values of K. There are two branches of excitation spec- 
trum: one is a gapless and linear mode at low momenta, 
denoted by (jj^{q), and the other is a gapful mode, de- 
noted by w^(q). Since w^(q) > until K exceeds a 
certain critical value, in the stability phase diagrams of 
Fig. [31 there is a region where the CSS state with finite 
superflow is stable. Like in the SF phase, we find two sce- 
narios regarding the instability of superflow in the CSS 
phase. First, sufficiently away from the boundary to the 
SS2 phase, Wj^(q) is pushed down with increasing K, and 
it reaches zero at K = K^^a, (see Figs.[4ljb) and (d)), sig- 
naling LI. In Fig. [21 -ft^cgg is plotted by the dotted blue 
lines separating the LI and stable CSS regions. When K 
is increased further, DI caused by phonons at low mo- 
menta sets in at 7f = K^gg as seen in Figs.^c), (c), and 
(f). K^gc, can be determined by the condition jfi^ — 0, 
which is reduced to 



1) 



{C^Vf{2cos{Kd) 
When C_y ^ J, we obtain 



J2(l + cos{Kd)f cos{Kd)(,23) 
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In Fig. [31 the right solid red lines represent K, 
is obvious from Eq. ([^3|) and Fig. [31 that K^gg is inde- 
pendent on n and monotonically increases with V/ J. It 
is also worth stressing that -R^^ss distinctively smaller 
than K^p. We attribute this reduction of the critical 
quasi-momenta to the difference between the energy band 
structures in the SF and CSS phases, which are respec- 
tively given by Eqs. ([T71) and ([^ . 

Secondly, near the boundary to the SS2 phase, DI 
caused by the excitations at q = (7r/d, 0) precedes the 
other instabilities and signals the transition to the SS2 
phase. The critical value of K for this DI is plotted by 
the dashed green line in Fig. [31^a). To complete the sta- 
bility phase diagram, we also carry out stability analyses 
for the SS2 state with superflow and locate the stable, LI, 
and DI regions. Notice that the LI region is almost in- 
visible in Fig. [31Ja) because the critical quasi-momentum 
for LI is nearly equal to that for DI. 

Finally, we discuss the feasibility of measuring experi- 
mentally the critical quasi-momenta. Experiments of ul- 
tracold dipolar Bose gases confined in a moving optical 
lattice will be performed in an additional parabolic trap. 
Since the density vanishes at the edge of the trapped gas, 
the critical quasi-momentum for LI is zero. Hence, for 
observing superflow, the temperature has to be so low 
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that the thermal component is invisible and LI cannot 
destabilize the system [23, IH, 113 ■ Another consequence 
of the parabolic trap is that the CSS phase inevitably 
coexists with other phases if it is present in the trap. For 
instance, when Uj = 0.4 at the trap center and V = 3.5 J, 
the CSS phase occupies the central region of the trap and 
is surrounded by the SF phase. Wc recall that the critical 
quasi-momcntum -ft^^gg for the DI in CSS is smaller than 
that in SF and independent on the density. This means 
that even in the trapped CSS phase coexisting with the 
SF phase, one can observe dissipationless flow by moving 
the optical lattice with a velocity smaller than K^a,a,/m 
and the breakdown of the superflow when KQgc^/m is ex- 
ceeded. Thus, given that the central density and V/J are 
precisely controllable in experiments, the stability phase 
diagrams of Fig. |3] can be investigated in current experi- 
mental setups. 

VI. SUMMARY 

In conclusion, we have studied stability of superflow of 
dipolar Bose gases in a moving optical lattice. Specifi- 
cally focusing on the superfiuid (SF) and checkerboard 
supcrsolid (CSS) phases, we have calculated the critical 
quasi-momenta for Landau and dynamical instabilities. 
Superflow in the CSS phases has been found to be sta- 
ble until the quasi-momentum exceeds the critical value, 
which is significantly smaller than that in the SF phase. 
In the CSS phase, we also found the dynamical insta- 
bility caused by roton-likc excitations that results in the 
transition to another type of supcrsolid. We emphasize 
that measuring the critical quasi-momenta will be a di- 



rect signature of superfluidity of the SS phases. 

Let us make brief comments on two recent experi- 
ments [13, [in that have explored supcrsolid phases in 
different contexts from dipolar bosons in optical lattices. 
One of them has studied spin textures in spin-1 Bose- 
Einstein condensates (BEC) of ^^Rb atoms with the fer- 
romagnetic contact interaction and the dipole-dipole in- 
teraction [soj . It was found that the spinor condensate 
has spatial magnetic order. However, this experiment 
is not convincing enough to proclaim the discovery of a 
supcrsolid phase not only because the magnetic order is 
short-ranged, but also because the equilibration time of 
the system is so long that one can not judge whether the 
magnetically ordered state is really an equilibrium state. 
The other experiment of Ref. [3l| has studied the system 
of a BEC coupled with an optical cavity and observed the 
formation of checkerboard density wave order in the BEC 
associated with the superradiant phase transition. Since 
the superfluidity of this possible supcrsolid phase has not 
been confirmed yet, it will be important to investigate 
the critical velocity in this system both theoretically and 
experimentally. 

Note added: After the submission of the present paper, 
there appeared relevant work by Kunimi et ai, which 
studies the critical velocit y of a supcrsolid phase through 
a single barrier potential |45| . 
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